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Abstract. To date, the quasar 1556+3517 is the only ra- 
dio loud Broad Absorption Line QSO (BALQSO; Becker 
et al. 1997). This prompted narrow band filter imaging 
observations in the range 4-15 fi with the ISO satellite. 
The source is clearly detected in all filters and appears 
point-like and isolated at the resolution of ISO. The over- 

\ all Spectral Energy Distribution (SED) in erg s _1 peaks 
at a rest wavelength > 6 fi and is reddened by 1.6 visual 
magnitudes. Correction for reddening brings 1556+3517 
within 1.3 sigma from the logR = 1 dividing line between 
radio-loud and radio-quiet sources. The mid-IR luminos- 
ity integrated over the 1.8-6 /i rest wavelength range is 
- 2 1 46 ergs- 1 (H = 75 kms" 1 Mpc" 1 , q = 0), which 

■ requires at least 4 Mq of dust and 800 Mq of associated 
gas. It is unlikely that such a large mass stems from the 
BAL wind itself. 
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1. Introduction 

Broad Absorption Lines (BAL) quasars constitute about 
~ 10 % of the overall quasar population. Observationally, 
they are characterized by very broad (0.1 c at zero inten- 
sity) and blue shifted absorption lines mostly from highly 
ionized species such as NVA1240, CIVA1550, .... Approx- 
imately 10 % of the BALQSO display in addition absorp- 
tion lines from low ionisation species such as Sill, Mgll 
and A1III (Weymann et al. 1991). These are the so-called 
low ionisation BALQSO, or Lo-BAL. BALQSO in general 
are also heavily absorbed in the X-rays (Green & Mathur 
1996) though it is not clear that the hot and highly ion- 
ized gas needed to produce the ionization edges seen in 
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the keV range is the same as that which produces the UV 
and optical absorption troughs. 

The BALQSO phenomenon is poorly understood the- 
oretically. In particular, its relation to the general QSO 
phenomenon is unclear. There are basically two classes 
of models: in the first class, BALQSO represent a transi- 
tory "cocoon" phase in an evolutionary process whereby 
a temporary excess of accretion fuel is expelled by radia- 
tion pressure from the central engine through a wind (Voit, 
Weymann & Korista 1993). This model predicts that even- 
tually BALQSOs will evolve into "normal" quasars once 
the excess gas has been removed. In the second class of 
models, BALQSO are normal quasars but viewed at a par- 
ticular inclination such that our line-of-sight intercept a 
wind of material blown-off the surface of the accretion disk 
or the top of the molecular torus (e.g. Murray et al. 1995). 
Spectropolarimetric observations tend to support this sec- 
ond hypothesis in that they often reveal significantly po- 
larized emission which is best explained by the scattering 
of a partially obscured but otherwise normal quasar spec- 
trum (e.g. Wills et al. 1992, Cohen et al. 1995, Hines & 
Wills 1995, Goodrich & Miller 1995 ...). The fact that 
the emission lines and continuum properties of BALQSOs 
are indistinguishable from that of non BAL QSOs (Wey- 
mann et al. 1991) also lend credit to the viewing direction 
hypothesis. 

An intriguing and important property of BALQSOs is 
that, up to now, there were no known radio-loud broad 
absorption line quasars (Stocke et al. 1992). The criterion 
commonly used for radio- loudness is logR > 1, where R is 
the K-corrected ratio of 5 GHz to 4400 A (B Band) flux 
(Stocke et al. 1992). There is some intriguing suggestion 
that BALQSOs show an excess of radio emission over non- 
BALQSOs such that there is an apparent accumulation 
of BALQSOs close-to but slightly below the logR = 1 
dividing line (Francis, Hooper and Impey 1993). 

The quasar 1556+3517 (z = 1.48) was recently dis- 
covered in the course of the FIRST VLA radio survey 
(Becker White and Helfand 1995). 1556+3517 belongs to 
the Lo-BAL class, and, with log R = 3.1, unambiguously 
qualifies as the first and so far only radio-loud BALQSO 
(Becker et al. 1997). Because it holds the promise of 
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shedding light on the radio- loudness/BALQSO dichotomy 
and, by extension, on the difference between radio-loud 
and radio-quiet quasars, I initiated a program to observe 
1556+3517 in the mid-IR with the ISO satellite (Kessler et 
al. 1996). One particular aim was to obtain a better esti- 
mate of the amount of dust on the line of sight to the QSO 
and check that 1556+3517 remains radio-loud even after 
reddening correction. Another was to verify the prediction 
that Lo-BAL QSO in general are heavily obscured. More 
generally, it is important to measure the spectral energy 
distribution (SED) of BALQSO over as wide a frequency 
range as possible in order to compare it to that of non- 
BAL QSO in the hope of finding differences that hold clues 
to the origin of the BAL phenomenon. 

2. Observations and data reduction 

The quasar 1556+3517 was observed with the 32 x 32 px 
long wavelength (LW) array of the ISOCAM instrument 
(Cesarsky et al. 1996) on board the ISO satellite (Kessler 
et al. 1996). The observation was executed nominally with 
the C01 observing mode (AOT) on 14 March 1997 (revo- 
lution 484), starting at 06:17:10 (UT) and lasting 46 min- 
utes and 26 second. It consisted of a series of 3 x 3 raster 
maps centered on the optical position of the QSO (R.A. = 
15h 56m 33.8s, Dec = 35° 17' 58", epoch = J2000), with 
one map per filter. The magnification selected was 3" per 
pixel and the distance between two raster points was 6". 
The axis of the rasters were oriented parallel to the celes- 
tial North and East. The number of exposures per raster 
dwell point was 19, with a unit integration time of 2.1 
s per exposure and an amplifier gain of 2. The 6 filters 
used were, in sequential order: LW3 [12-18 /i], LW2 [5- 
8.5 n], LW9 [14-16 /j], LW8 [10.7-12 fj], LW7 [8.5-10.7 fj] 
and LW1 [4-5 /z], where the numbers in brackets refer to 
the wavelength range covered by each of the filter. This se- 
quence corresponds to decreasing pixel illuminations such 
that it minimizes the detector stabilization time. 

The six monochromatic maps were processed with the 
CAM Interactive Analysis (CIA) Software 1 starting from 
the Edited Raw Data (ERD). The first few (~ 4) expo- 
sures of each frame were discarded so as to retain only 
those data for which the detector throughput had reached 
90 % of its stabilisation level. The remainder of the pro- 
cessing was standard and involved the usual steps, i.e. 
dark current subtraction, flat-fielding, removal of particle 
hits ( "de-glitching" ) and flux calibration. The source is 
clearly detected in all six maps at a very high significance 
level (e.g. 63 sigma for the LW3 map). The 12-18 /j map 
is shown in Fig.l as an example. 

The monochromatic intensities were computed by in- 
tegration of the flux in a 6 px radius (18") circle centered 

1 CIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOCAM Consortium led by the ISOCAM PI, C. 
Cesarsky, Direction des Sciences de la Matiere, C.E.A., France. 




Fig. 1. The 12-18 /i map of the field around the quasar 
1556+3517; the size of the map is 75 x 84 " and the scale 

3" pic^- 
on the source and subtraction of the normalized back- 
ground integrated in an annulus of 8 and 13 px inter- 
nal and external external radii, respectively. The results 
are displayed in table 1. The errors quoted derive mainly 
from the uncertainty in the calibration and background, 
the latter being induced by residual flat-fielding errors. 
The full widths at half maximum (FWHM) of the source 
images were computed by fitting a two dimensional gaus- 
sian to the flux spatial distribution. The results are listed 
in column 5 of table 1. For comparison, the theoretical 
value of the full width at 41 % of the maximum intensity 
of a diffraction limited telescope of diameter D is p, i.e. 
0.343 A (/i) in the case of ISO. This yields e.g. 5.0" for 
the LW3 image, consistent with the value in table 1. As 
a sanity check, the LW3 images of point sources obtained 
with the same pixel magnification were retrieved from the 
set of ISO calibration files and analysed in the same fash- 
ion, yielding FHWM =4.4 ", in good agreement with the 
FWHM of 1556+3517. For all six filters, the source width 
is consistent with that of a point source. Neither of the 
six maps show extended structures nor close-by objects 
at least down to the detection limits quoted in column 6 
of table 1. These limits were obtained by computing the 
standard deviation of the background in areas of the map 
well outside the image of the source. 



3. Discussion and conclusions 

The overall rest frame Spectral Energy Distribution (SED) 
of 1556+3517, from 1.5 GHz to the B band, is plotted in 
fig. 2, where the optical and radio data are taken from 
Becker et al. (1997). All observed fluxes have been con- 
verted to rest frame luminosities vL ve (in ergs -1 ) assum- 
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Table 1. Monochromatic mid IR fluxes of 1556+3517 



LW 


AA 


^rest 


Flux 


FWHM 


Limit 






(aO 


(mJy) 


(") 


(/iJy/n) 


1 


4.00-5.00 


1.81 


1.5 + 0.4 


3.5 + 1.0 


14 


2 


5.00-8.50 


2.75 


4.7 + 0.5 


3.0 + 0.4 


8 


7 


8.50-10.7 


3.90 


5.7 + 0.3 


2.9 + 0.5 


11 


8 


10.7-12.0 


4.56 


9.8 + 0.5 


3.2 + 0.6 


20 


3 


12.0-18.0 


5.85 


12.7 + 0.5 


4.3 + 0.5 


8 


9 


14.0-16.0 


6.01 


8.3 + 1.5 


3.9 + 1.0 
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ing H = 75 km s 1 Mpc 1 , q = and using the standard 
equations (e.g. Weedman 1986): 



L ve = Aird 2 {1 + z) f vo 



d- 



cz(l + z/2) 
H (1 + z) 



(1) 
(2) 



where L„ e is the monochromatic luminosity per unit 
frequency in the rest frame of the QSO, i vo is the 
monochromatic flux per unit frequency in the observer's 
frame, d is the distance to the QSO and the other symbols 
have their usual meanings. 

The SED of the quasar IRAS 14026+4341 is also 
shown in fig. 2 for comparison, the data being taken from 
Low et al. (1989). The comparison is illustrative since 
IRAS 14026+4341 is an IRAS discovered BALQSO with 
a large IR-to-optical luminosity ratio (Low et al. 1989). It 
is immediately apparent that, despite similar ultraviolet 
luminosities, 1556+3517 is almost an order of magnitude 
more luminous in the near to mid-IR than 14026+4341. 

This presumably indicates the presence of a large 
amount of warm dust in 1556+3517 which causes a de- 
pression of the optical-U V continuum by reddening and 
an enhancement of the near and mid-IR flux because of 
grain thermal emission. 

The 1556+3517 spectral index (F„ cx v a ) from 
the B band (X rest = 0.1774 /z) to the LW1 filter 
(\ rest = 1.81 n) is a = -1.99 ± 0.25. Using the R band 
(\ rest = 0.2826 n) flux instead yields a = -1.46 ± 0.31. 
Both indices are significantly steeper than those of "nor- 
mal" non-BAL quasars. Note that only the hottest dust 
grains at ~ 1600 K yield a significant contribution to the 
1.81 fj, flux. Since such a temperature is higher or equal to 
the grain sublimation temperature, relatively little dust 
thermal emission is expected at this wavelength. Hence, 
to a first approximation, the steep UV-to-near-IR spec- 
tral index of 1556+3517 is mostly a consequence of dust 
extinction. 

There are two reasons why the a = —1.99 + 0.25 index 
value should be preferred: 

1. The 1774.2 A band lies in a region relatively free of 
any strong emission or absorption lines. This is not 
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Fig. 2. The overall spectral energy distribution (SED) of 
1556+3517 in the QSO rest frame from radio to ultraviolet. 
The SED of the IR selected BALQSO IRAS 14026+4341 is 
also plotted (circled symbols connected by a dotted line) for 
comparison. 



the case of the 2823 A band which is contaminated by 
the MgUA 2798 emission and absorption. 
2. More importantly, the 2823 A band is definitely en- 
hanced by Balmer continuum (BaC) emission and 
blended low-contrast Fell emission lines which create 
a pseudo-continuum in quasar spectra often referred 
to as the "small-bump" . 

The extinction internal to 1556+3517 was therefore 
computed under the following assumptions: 

1. The "intrinsic", i.e. un-reddened, optical-UV spec- 
tral index of 1556+3517 is -0.5. Within the uncer- 
tainty, this is equal to the mean UV spectral index 
a = —0.46 ± 0.05 of the sample of 33 high red- 
shift (z> 1) QSOs of O'Brien, Gondhalekar and Wil- 
son (1988) as well as to the mean optical spectral in- 
dex a = -0.5 ±0.1 of the large 227 QSO sample of 
Cheney and Rowan-Robinson (1981) or that obtained 
by Neugebauer et al. (1987) from their sample of 104 
quasars (a = -0.4 + 0.04). 
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2. The reddening law is as parametrized by Seaton (1979) 
in the optical-ultraviolet range and as tabulated by 
Ricke and Lebofsky (1985) from 1 to 13 \i. 

Under these conditions, the amount of visual extinc- 
tion internal to 1556+3517 required to steepen the power- 
law index from -0.5 to a = —1.99 is given by: 

(a + 0.5) Zn^) 

A - ' v A 2 ' (a) 

Jn(10°- 4 («*i-Kx a ) (6) 

where [Ai — A2] is the (rest) wavelength range over 
which the spectral index is measured (0.177-1.810 /i in 
the present case) and the function R(X) = -j^ is the red- 
dening law at wavelength A. This yields a visual extinc- 
tion A v = 1.64 ± 0.27 magnitudes which implies that the 
continuum is reddened by 4.01 magnitudes at a rest wave- 
length of 1774 A but only 0.256 magnitudes at 1.81 /U. 
Note that the foreground galactic extinction in the di- 
rection of 1556+3517 (In = 56.48°, b n = 49.95°) is only 
0.05 V magnitudes and can be neglected. 

After correction for extinction, the radio-to-optical 
flux ratio becomes logR = 1.50 ± 0.23 if the B band data 
are used to derive the flux at A res t = 4400 A (K correc- 
tion) and logR = 1.15 + 0.25 if one uses the R flux in- 
stead. The two estimates differ at the 1.1 sigma level only. 
I therefore adopt the mean of the two log R ~ 1.32 + 0.25. 
Such a ratio is less than 1 .3 sigma away from the log R = 1 
dividing line between radio-loud and radio-quiet quasars. 
One therefore concludes that the case for 1556+3517 
being radio-loud is marginal, at best. It is also worth 
pointing-out that de-reddening brings 1556+3517 close to 
the log R = 1 line where data suggests that there is an ap- 
parent excess of BALQSOs over non-BALQSOs (Francis, 
Hooper and Impey 1993). 

The near to mid-IR luminosity Lmir of 1556+3517 was 
computed by integrating all the de-reddened rest-frame 
emission over the range A rcs t = 1.81 /i to A ros t = 6.01 /1 
after subtraction of an -0.5 index power-law matching 
the de-reddened B band and LW1 band luminosities: 
Lmir = 1.9 10 46 ergs -1 . 

One can use Lmir to estimate the mass of dust re- 
sponsible for near to mid-IR emission in 1556+3517. The 
mean rest wavelength of the MIR emission corresponds to 
a black-body temperature of ~ 740 K. Under the assump- 
tions that the dust grains have a mean size of 0.1 and 
radiate like black-bodies and following Rudy and Puetter 
(1982), I obtain a dust mass M d ~ 4Mq. This should in 
fact be seen as a lower limit since actual grains are likely to 
emit less efficiently than black-body thereby requiring an 
even larger mass to produce the observed MIR luminosity. 
For a normal galactic dust-to-gas mass ratio of 200, this 
implies a total mass of gas of about 800 Mq . 

The color temperature of the dust is higher than that 
of typical HII regions which suggests that it is closely as- 
sociated to the active nucleus. It is therefore legitimate to 



ask whether the dusty gas responsible for the MIR emis- 
sion could be the BAL wind itself. In the disk wind model 
of Murray et al. (1995), the BAL region covering factor is 
about 10 %, its column density Nh ~ 10 23,5 cm -2 and its 
inner radius is of the order of 10 16 cm. This yields a BAL 
mass Mbal — 3 10 -2 Mq, more than 4 orders of magni- 
tude smaller than the above estimate. Of course, there is 
more mass available as the wind will accumulate material 
over the years. For an axially symmetrical wind, the mass 
flow rate is given by: 

M = 3.33 10- 2 {—)N 23 R 16 v 8 (M© yr- 1 ) (4) 

where, ^- is the solid angle sustained by the wind as seen 
from the central source, N23 its column density in units 
of 10 23 cm -2 , Ri6 its radius in units of 10 16 cm, and vg its 
velocity in units of 10 8 cms -1 . Assuming parameters value 
as in Murray et al's model yields M ~ 3.3 10~ 3 Mq yr -1 . 
This implies that the BAL phenomenon should have lasted 
at least ~ 250,000 years to accumulate the MIR emitting 
gas. While such a value in itself is not unreasonable, there 
are problems with this scenario: 

— On the one hand, the dust must somehow "see" the 
central UV source in order to retain a temperature of 
~ 700 K for such a long time. However, given the mass 
involved and the geometry of the wind, the opacity will 
be extremely high and the bulk of the dust and gas will 
be (self-) shielded from the central source. 

— This gas would presumably be at least partially 
ionized. Since the total column would far exceed 
10 24 cm -2 , the gas would be opaque to free-free absorp- 
tion and the continuum totally suppressed, contrary to 
what is observed. 

It therefore seems unlikely that the observed MIR 
emission originates from the wind itself. A more plausi- 
ble origin would be e.g., the molecular torus, the Narrow 
Line Region, or a recent burst of star formation. 
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